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We discuss a hybrid quantum system where a dielectric membrane situated inside an optical cavity 
is coupled to a distant atomic ensemble trapped in an optical lattice. The coupling is mediated by 
the exchange of sideband photons of the lattice laser, and is enhanced by the cavity finesse as well as 
the square root of the number of atoms. In addition to observing coherent dynamics between the two 
systems, one can also switch on a tailored dissipation by laser cooling the atoms, thereby allowing 
for sympathetic cooling of the membrane. The resulting cooling scheme does not require resolved 
sideband conditions for the cavity, which relaxes a constraint present in standard optomechanical 
cavity cooling. We present a quantum mechanical treatment of this modular open system which 
takes into account the dominant imperfections, and identify optimal operation points for both 
coherent dynamics and sympathetic cooling. In particular, we find that ground state cooling of a 
cryogenically pre-cooled membrane is possible for realistic parameters. 
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I. INTRODUCTION 

Micromechanical resonators coupled to cold atoms 
have recently attracted much interest as a promising class 
of hybrid quantum systems ■1-3]. On the one hand, ad- 
vances in nanoscale fabrication techniques have led to the 
development of high-quality micromechanical resonators 
in various forms, driven in particular by the interest in 
preparing quantum states of macroscopic devices [1] . An 
attractive feature of these resonators is that they may be 
functionalized by attaching additional elements, which 
enables optical, magnetic, or charge-based coupling to 
other systems [3]. Many such devices are optimized for 
the interaction with light in cavity-optomechanical se- 
tups [MS]: which, e.g., enable optomechanical laser cool- 
ing [9HT8] of the mechanical resonator to its ground state 
[TM21] . In addition to being useful for applications in 
metrology [25], the control of these devices on a single 
quantum level will also enable addressing questions in 
quantum foundations [4, 23| 124] and applications related 
to quantum information [25H34] . On the other hand, cold 
atomic gases are systems over which an unprecedented 
level of quantum control has been reached, involving both 
internal and external (motional) degrees of freedom |35j . 
Here, a comprehensive toolbox for cooling, state prepa- 
ration and readout is available, profiting from ongoing ef- 
forts in metrology [36, 37 and quantum simulation 38J. 
These efforts have also enabled an atomic realization of 
optomechanics [55H4T] . 




FIG. 1. (a) Micromechanical membrane in a cavity coupled 
to a distant atomic ensemble by means of a laser (red, fre- 
quency The motion of the membrane shifts the opti- 
cal lattice dipole trap of the atoms, while the motion of the 
atoms changes the radiation pressure on the membrane. In- 
set: The atoms are modeled as two-level systems, which are 
far off-resonant with respect to the laser, (b) Refractive index 
profile used in the one-dimensional model (Eq. (pi). 



In light of these developments, a combination of me- 
chanical and atomic systems promises new ways to con- 
trol and manipulate the mechanical resonator via the 
atoms, but may also lead to novel ways of realizing in- 
teresting atomic dynamics. Along these lines, direct 
[2j[3] and optical [42H48] coupling schemes between atoms 
and micromechanical systems have been proposed, and 
first experiments have been realized [4"5ll51| . Many of 
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these setups pose great experimental challenges as they 
require the combination of a functionalized mechanical 
resonator with atomic trapping and manipulation tech- 
niques within the same vacuum chamber or cryostat and, 
possibly, within the same optical cavity. 

A modular setup that circumvents the need for inte- 
gration of atoms and mechanics within the same vac- 
uum chamber has recently been proposed :48J and later 
realized [5T] by some of the authors. The experiment 
demonstrates a light-mediated coherent coupling of the 
vibrations of a micromechanical membrane to the center- 
of-mass motion of an atomic ensemble, which is kept at 
a distance (1 m in the experiment [51] ) in a separate vac- 
uum chamber. Motivated by these experimental develop- 
ments, we analyze in this work a more sophisticated setup 
as displayed in Fig.JIja). Here, the membrane is enclosed 
in an optical cavity, as in the standard 'membrane-in-the- 
middle' configuration |16j . while the atoms are trapped 
outside of the cavity in the optical lattice potential pro- 
vided by the light field reflected off the cavity. Adding 
the cavity preserves the modularity, while the coherent 
coupling between the two systems is enhanced by the cav- 
ity finesse in addition to a collective enhancement by the 
square root of the number of atoms. As we will show, 
the coupling strength can be leveraged so as to exceed 
the various decoherence rates. First of all, this coherent 
coupling enables a resonant exchange of excitations be- 
tween the two systems. Furthermore, one can switch on 
laser cooling of the atomic motion [52] , thereby effectively 
tailoring the dissipation of the membrane. This leads to 
a sympathetic cooling scheme for the membrane which 
is analogous to standard optomechanical cavity cooling 
[H [10], with the role of the cavity played by the atoms. 
However, in addition to being tunable, our scheme does 
not require the cavity linewidth to be much less than the 
mechanical frequency and thus relaxes the experimental 
requirements on the cavity. 

The aim of this work is to provide a rigorous, quantum 
mechanical treatment of the proposed setup. In addition 
to deriving the exact form of the coherent coupling, this 
also yields a consistent description of quantum noise in 
the problem such that we are able to describe the dynam- 
ics of the system on the level of a few quanta. In par- 
ticular, we analyze the possibility of observing quantum 
coherent excitation exchange between the two systems, 
and address the question whether the membrane can be 
cooled to its quantum mechanical ground state by laser 
cooling the atoms. Our results allow us to identify opti- 
mal experimental parameters for both scenarios. 

The remainder of this paper is structured as follows: 
In Sec.|H] we first give a simple semi-classical estimate 
for the coupling, and then present a full quantum model 
of the problem from which we obtain a linearized in- 
teraction of atoms and membrane with the light field. 
Section III then deals with eliminating the light field to 



discussed in Sec.|IV[ where we analyze the scenarios of 



obtain an effective description of atoms and membrane 
alone, which includes the dominant sources of decoher- 
ence. The resulting effective dynamics is subsequently 



coherent dynamics (Sec. IV A) and sympathetic cooling 
of the membrane (Sec jlVB ). Finally, Sec. [V] comments 
on the experimental realization of our proposal, and con- 
cluding remarks can be found in Sec.|VI| 



II. MODEL 

We consider a system as shown in Fig.[IJa), where a 
mechanical resonator residing in a single-sided cavity of 
finesse J- is coupled to the motion of a distant atomic 
ensemble outside the cavity. Before discussing our full 
quantum model below, we present a brief estimate of the 
expected coupling based on a quasi-static picture: the 
membrane is taken to be a standard SiN membrane with 
a resonance at frequency w m and effective mass M , such 
that its zero point motion is given by l m — yJh/Mbj m . 
Its coupling to the atoms is mediated by a laser beam of 
wavelength Al , which is reflected off the cavity and forms 
an optical lattice for the atoms by means of the dipole 
force [53]. We assume the lattice to be deep enough to 
treat its wells as independent harmonic traps of frequency 
uj at ~ u) m . If the laser is resonant with the cavity, then a 
displacement of the membrane by an amplitude l m cor- 
responding to its zero point motion leads to a phase-shift 
8<j) ~ J-lm/^L of the reflected light. For a single atom, 
the resulting shift of the standing wave pattern outside 
the cavity would lead to a force oc muj^Scjj/k^, where 
m is the mass of one atom and k^ the wavenumber of 
the laser. However, since all atoms are affected equally 
the membrane couples to the center-of-mass motion of 
the atoms, which, on a single quantum level, leads to a 
collectively enhanced coupling rate of the form [33 1ST] 
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Here, N is the number of atoms and l at — -\/fi./mw a t is 
the atomic zero point motion. This estimate shows that 
the small mass ratio m/M can be compensated by both 
a large atom number and a high finesse. In reverse, if the 
atoms depart from their equilibrium positions, then their 
restoring force originates from redistribution of photons 
between the right- and left-propagating components of 
the standing wave. As a result, the intensity inside the 
cavity changes and thus the radiation pressure on the 
membrane such that it experiences a force. The quantum 
model presented in the following will confirm the above 
estimate. 



A. Quantum Description 

We describe our setup by a one-dimensional (ID) 
model with Hamiltonian 



H — H m + H at + -ffficid + H m .{ + H at .f + iJdcc • 



(2) 
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Here, H m , H at and -fffieid describe the free evolution of 
the mechanics, the external atomic degrees of freedom, 
and the field modes, respectively. Further, H m _ { repre- 
sents the coupling between the mechanical oscillator and 
the field modes, and H at -f the interaction between atoms 
and field modes. The last term in Eq. ([2| summarizes 
decoherence effects, which will be introduced explicitly 
in SecHEE] 

The first term in Eq. ([2| describes the harmonic vibra- 
tions of the mechanical resonator and reads 



2M 



hu>„ 



(3) 



where z m , p m and a m are the position, momentum and 
annihilation operators, respectively. The second term in 
Eq. ([2]) is given by H at = Ylj=iPj an< ^ represents 
the kinetic energy of the atoms, where the pj are the mo- 
menta of the individual atoms. The free field Hamilto- 
nian reads Hg e id = / dw fkv b^b^ , where the mode opera- 
tors obey the commutation relation \b u ,b^,] = S(uj — oj'), 
and all frequency integrals are understood to run over a 
bandwidth 29 around the laser frequency ujl introduced 
below. In the spirit of a 'modes of the universe' approach 
we take care of the cavity by explicitly including it in the 
mode functions u u (z), which are associated with the b u 
and defined by writing the positive frequency parts of the 
fields as [53] 

E (+) (z) = i J du>S u u u {z)b u , (4) 
P(z) = ^d W U(*. (5) 



Here, — \J huj / ire^cA, u u (z) = d z u u {z) jik and A is 
the cross sectional area of the modes. In the following we 
discuss the particular case of a 'membrane-in-the-middle' 
setup [TB] based on a single-sided cavity, as shown in 
Fig.jlja). However, other arrangements can be treated 
analogously by using the appropriate mode functions. 
This is shown explicitly in Appendix [A] for a setup with 
a movable end mirror instead of the membrane. 

We model the membrane and the semi-transparent end 
mirror of the cavity by dielectric slabs at positions i and 
L, respectively [THIES]- The wave equation for the mode 
functions thus has to be solved for the refractive index 
profile 



n(z) 



if ze [1,1 + d m ] 
if z e [L, L + d] , 
otherwise 



(6) 



which is illustrated in Fig.Qb), and with the boundary 
condition of a vanishing field at z — 0. The latter cor- 
responds to the limiting case of an ideal mirror at posi- 
tion 2 = 0. In the above equation, thickness (d, d m ) and 
refractive index (n, n m ) of each slab determine the asso- 
ciated amplitude reflection and transmission coefficients, 



which read 



t 



(n 2 — 1) sin(kdn) 



2incos(kdn) + (1 + n 2 ) sm(kdn) 
2in 

2incos(kdn) + (1 + n 2 ) sm(kdn) 



(7) 
(8) 



for the end mirror and analogous expressions hold for the 
coefficients r m and t m of the membrane. As a result, the 
properly normalized mode functions read 



u u (z) = 




-,—ik 



-ikz 



if < z < i 

if £ + d m < z < L 
if z > L + d 



where we have defined the coefficients 
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(9) 

(10) 
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Clearly, ip u is the phase by which the standing wave pat- 
tern outside the cavity is shifted. In the limit of an almost 
transparent membrane and a high reflectivity end mirror, 
given by t m — > 0, |t| — > 1, the factor T u represents the re- 
sponse of an empty single-sided Fabry-Perot cavity with 
well-defined resonance frequencies cj„ = vkc/L 1 where 
v = 1,2,.... Around the latter, |T W | 2 can be approxi- 
mated by a Lorentzian, i.e. 



2T 



IT K 2 + (UJ — UJ U ) 2 



(12) 



Here, k = c(l— |r|)/2i|r| 1 / 2 is half of the cavity linewidth 
and JF = ttc/2kL is the finesse, which we assume to be 
limited by the semi-transparent end mirror. For small 
displacements z m the radiation pressure interaction be- 
tween membrane and field in Eq. ^ is given by [481 157j 

e K - l)A 



E(-\l + d nl )E(+\l + d m ) 
-E^(£)E^(£)]z m , (13) 



which we will evaluate below. 

Finally, H at -t in Eq. ^ describes the dipole interaction 
of the atoms with the light field. We model the atoms as 
two-level systems of transition frequency w eg , from which 
the frequency wl of the trapping laser is far detuned by 
6 = WL-w eg (see inset in Fig.[lja)). We further anticipate 
that all relevant field modes will lie within a bandwidth 
29 <C 1 8 1 around the laser frequency, such that we can 
eliminate the upper level of the atoms and describe the 
atom-field interaction by an AC-Stark shift for atoms in 
the ground state 



at-f 



HS 



N 



(14) 
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where /i is the atomic dipole moment and the Zj are the 
positions of the atoms. 



with the phase difference Atp u — ip u — ip^ L and the 
membrane-field coupling constant 



B. Linearization Around Laser Drive 

We now introduce the laser displayed in Fig.JIJa), 
which provides the optical lattice for the atoms and medi- 
ates the coupling between the atoms and the mechanical 
resonator. To do so, we carry out the replacement for 
the field mode operators 



(15) 



where the amplitude a is related to the running wave 
power P — hw-^a 2 / 2tt of the incident laser. In the fol- 
lowing we assume a 3> 1, which allows us to linearize 
the interaction Hamiltonians i? m -f and i? a t-f around the 
laser by only keeping terms of order a and a 2 . 



1. Membrane-Field Interaction 

We start with the membrane-field interaction H m .{ by 
inserting the field expansion in Eq. Q together with the 
mode functions given in Eq. ^ into Eq. (13 1, and then 
apply the above replacement to the field operators. The 
contribution oc a 2 yields a constant force on the mem- 
brane that can be accounted for by a redefinition of 
its equilibrium position. The contribution linear in a 
provides the relevant interaction between mechanics and 
light field. For convenience we introduce A; jW = u u {l) 
and Xr.u, = u u {t + d m ) for the mode functions such that 
in a rotating frame with respect to the laser frequency 
the linearized membrane-field interaction H m _{ is given 
by the contribution linear in a: 



g m . f = co(< : 1) - 4a /,L- 



{A* A r>a , L — Af w A;, WL }& w 



At this point we need to make assumptions on the 
frequency-dependence of the mode functions u u (z) in or- 
der to proceed. We focus on the case of small membrane 
reflection |t m |, such that the cavity response is not al- 
tered significantly by the presence of the membrane and 
assume that the laser is resonant with a cavity mode. 
As motivated in Sec. |III| below, we further assume a 'bad 
cavity' regime, where the cavity half-linewidth k is much 
larger than the separation of the relevant sidebands from 
the laser frequency wl (see Fig.|2b. In this case we may 



approximate \T U 
by using S u 



|T WL |, such that u u (z) 



Xz), and 



£ Uh we obtain 



Hjn-i « hg n 



(b u e iA ^+ble- iA ^)(a m + al) 



(17) 



9m = - l)Aa£ 2 L U\\ r ^\ 2 - |A^ L | 2 ) . (18) 



In the limit k-^d m —> of a thin membrane and small but 
constant reflection |r m | « k-^d^n^ — l)/2 we can write 
the membrane-field coupling constant as 



ffm = ^n |rm|sin(2W)|T |2 
V 7r 

ak h l m 2T 
= — r~ Fm — , 



(19) 
(20) 



where we assumed sm(2ki J £) = 1 in the second line. This 
corresponds to placing the membrane on the slope of the 
intra-cavity intensity, where the coupling is maximal. In 
contrast, at points of zero intensity the coupling vanishes 
and higher orders in z m would have to be considered [16] . 
However, these points will not be relevant in this work. 



Note that the coupling in Eq. ( 20 ) is enhanced by the 



finesse and that a similar calculation for a cavity with 
movable end mirror yields exactly the same interaction in 
Eq. (17), but with g m — > g' m = 2ak\,l ra iF /tt 3 / 2 , as detailed 
in Appendix |A"| 



2. Atoms-Field Interaction 



The next step is the expansion of the AC-Stark shift 
Hamiltonian H at _{ in powers of a. We insert Eq. Q 
with the mode functions of Eq. |9| into Eq. ( 14 ) and ap- 
ply the replacement Eq. ( 15 ). The resulting contribution 
proportional to a 2 provides an optical lattice trap that 



is seen by each atom and that is described by the po- 
tential V(zj) — Vo sm 2 (k^Zj + (fu h ) with lattice depth 
Vo = [J, 2 £ 2 \a\ 2 /h5. The minima Ej for blue detuning 



(16) ^ 0) are determined by the condition k\^Zj-\-Lp u 



where v = 1,2,.... Combined with the kinetic energy 
of the atoms in iJ at a harmonic approximation around 
the trap minima defines the atomic trap frequency w at 
via muj 2 t /2 = Vok 2 ^. The contribution linear in a yields 
the interaction between light field and atomic ensemble. 
By applying a Lamb-Dicke expansion around the equilib- 
rium positions of the atoms Zj, we obtain a linear inter- 
action between the field modes and the fluctuations 
around the equilibrium position of each atom given by 
z j = lat(a) + o>j)/V2, where [ctj, at] = 1. In carrying out 



the remaining sum in Eq. ( 14 1 we replace the zj by the 



center-of-mass position z of the whole ensemble, which 
is a good approximation if w a tL at /c <C 1, where L at is 
the linear extension of the ensemble. For typical experi- 
mental parameters this assumption is fulfilled in general, 
see Sec.|V] As a result, the light field only couples to the 
atomic center-of-mass mode a a t = Si a j/V^V according 
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to 



Sat , .9m 



duj 



(b u e iA ^+bie- iA ^)(a at + al) 



x ship- k L )z + A(p u ], (21) 

where the atom-field coupling 

w at 



ffat 



h:\fN 



(22) 



2ak L l at 

is enhanced by the square-root of the number of atoms. 

3. Final Linearized Hamiltonian 

In summary, the final linearized Hamiltonian contains 
interactions between the field modes and the mechanical 
and atomic center-of-mass motion. We neglect the con- 
stant contribution in Eq. ^ and together with Eq. (17) 
and Eq. (21 ) the linearized Hamiltonian in a frame rotat- 



ing at the laser frequency reads 



+ huj at a L a at + / dw SA U b\,b u 



(23) 



%at 



duj 
V2tt 



+bl){a at + al t )sm(A w z/c+Aipu). 



Here, A w = lj — cjl and we have absorbed the phase fac- 
tors e ±lAv ^ appearing in Eqs. ( 17) and (21 ) into the field 



operators . The effect of the cavity response T w enters 
both via the membrane-field coupling <? m , which depends 
on |T WL | , and the phase ip u = arg{T w } contained in 
Atpu,. 



III. EFFECTIVE DYNAMICS 

We are interested in deriving the effective equations of 
motion (EOMs) describing the coupling between atoms 
and mechanical oscillator, which are obtained by elimi- 
nating the light field in a Born-Markov approximation. 
Motivated by the experimental considerations presented 
below in Sec.|VJ we focus in the following on the 'bad 
cavity' regime, where the cavity half-linewidth k is much 
larger than the atomic and mechanical frequency, i.e. 
k w at , u; m as depicted in Fig. [2] Also assuming the 
weak-coupling conditions uj at ,uj m 3^5at> ffatffm we 



find that the resonant terms in Eq. ( |23| involve field op- 
erators at the sideband frequencies uj = wl ± w a t,m) and 
we choose w a t ~ to ensure resonant interactions. We 
further assume that the laser is resonant with a cavity 
mode, i.e. wl ~ u> v as shown in Fig.|2j which yields 
a high intra-cavity intensity, but also ensures that both 
sidebands enter the cavity equally well, i.e., that there 
is no spectral filtering. In the described regime the cav- 
ity adiabatically follows the mechanical oscillator and the 



\ < 2k 



FIG. 2. Light fields involved in mediating the effective cou- 
pling. The laser-enhanced interactions (coupling constants 
g m and g a t) correspond to a conversion between laser pho- 
tons wl (black) and sideband photons ojl ±w a t,m (blue/dark 
and red/light). The response profile of the cavity is indicated 
by the green dashed line and easily accommodates both side- 
bands in the 'bad cavity' regime discussed in the text. 



atoms and needs not to be treated as an independent de- 
gree of freedom. The opposite 'resolved sideband' limit, 
where k <C w m , u; a t , would require a different treatment 
and is not considered in this work. Finally note that 
the assumptions described here already entered in the 
derivation of Eq. ( 23 1 . 



A. Elimination of the Light Field 

We start by stating the complete hierarchy of 
timescales for our calculation 

5 > 9 > -,k > w at ,w m > g^,gl t ,g m g a t, (24) 

T 

where r = z/c is the propagation time between mechan- 
ics and atomic ensemble. The left-most inequality is 
needed to model the atom-light interaction by an AC- 
Stark shift as described in Sec. lII Al We assume a broad 
spectrum of field modes such that its bandwidth 26 is 
large compared to all frequencies associated with the sys- 
tem evolution. Below, this allows us to distinguish the 
temporal order of the light field's interactions with the 
atoms and the mechanical resonator. Motivated by ex- 
perimental parameters (see Sec.|V]), we assume that the 
retardation time between the two systems is short com- 
pared to the system's evolution, i.e. 1/t> u) a t,u) m . Fi- 
nally, we have the bad cavity assumption as well as the 
weak-coupling assumptions w at ,w m > g% t , g^, g ai g m as 
discussed in the the previous paragraph. 

To derive the effective EOMs for mechanical oscilla- 
tor and atoms we start with the Heisenberg EOMs for 
the quadrature operators obtained from the linearized 
Hamiltonian given in Eq. ( 23 1 (see Appendix [B] for de- 
tails). We then formally integrate the EOMs for the 
field quadratures. By inserting the resulting expressions 
into the Heisenberg equations for mechanics and atoms 
and applying the Born-Markov approximation, we ob- 
tain effective quantum Langevin equations (QLEs) with 



G 



interactions retarded by a time f = r + l//c. The ad- 
ditional contribution 1/k represents the retardation due 
to the cavity and follows from the Taylor expansion of 
Atfu around the cavity resonance. From Eq. (24 1 we see 
that w a tT <C 1, so that we can neglect these retardations 
in the following, and thus obtain effective QLEs that are 
local in time: 



Pm 

x m 

Pat 



-uj m x m + gx at + \Jlnl S F m (t) 



— ^at^at 
WatPat ■ 



(25a) 

(25b) 
(25c) 
(25d) 



Here, F m (t) and 7^ lff are defined below, and g — 2g at g m 
is the effective coupling between the two systems. From 
the above equations we can read off the corresponding 
effective Hamiltonian 

H cB = t\h) Ia a\ a a m + fta> at <4 t a at - hgx m x at , (26) 

where the effective coherent coupling g contains the de- 
sired enhancements and is explicitly given by 



9 = 2g a t3n 



=^at^ — Fn 
tat. 



2F 



(27) 



The noise operator F m (t) appearing in Eq. (25a) is Her- 
mitian and represents the radiation pressure noise origi- 
nating from the intra-cavity light field, ft is characterized 
by the two-time correlation function 



(F m {t)F ID (t')) = 6e{t-t') ) 
and the associated diffusion rate 



(28) 



(29) 



appearing in Eq. (25a). The expression Sg(t — t') denotes 



a representation of the 5-function that is peaked on the 
timescale 1/9. As explained in Appendix [5] this momen- 
tum diffusion rate appears as a fundamental limitation in 
any optomechanical system, and it was recently observed 
in an experiment [58] . 



B. Decoherence and Imperfections 

In the following we will discuss the main decoherence 
channels of the mechanical oscillator and the atoms. In 
addition to the radiation pressure noise derived above, 
these include thermal decoherence due to finite support 
temperature and absorption of laser photons in the mem- 
brane, as well as light-induced momentum diffusion of the 
atoms. 

The mechanical resonator is coupled to a support of 
temperature Tq, which leads to a finite linewidth 7 m due 
to clamping losses, but also to thermal decoherence at 
rate y m N m s» k j^Q° , where ks is the Boltzmann constant 



and Q m = ui m /~/ m the mechanical quality factor. In addi- 
tion, the membrane is also heated by absorbing photons 
from the coupling laser. We model this heating effect by 
introducing an effective bath temperature T c r = Tq + 5T, 
where ST = is the temperature increase due to ab- 
sorbed laser power P a b s in the center of the membrane, 
representing a worst-case estimate |45) . The thermal link 
K t h connects the center of the membrane to the support 
and depends on the thermal conductivity and the thick- 
ness of the membrane, as well as the beam waist of the 
laser (see Appendix [C]). Further, the absorbed power for 
a membrane sitting on a slope of the intra-cavity inten- 
sity is given by F a bs = »mf^ where = I — - if n is 
the power absorption coefficient. The resulting effective 
bath occupation is thus given by 



N„ 



k B 



T 



ai4F_ 

Kth 7T 



P 



(30) 



as described in more detail in Appendix [C] 

Atoms in a 3D standing wave undergo diffusion pro- 
cesses. In our treatment these diffusion processes drop 
out as an artifact of the ID model, as already men- 
tioned in Appendix [B] and in Ref. [48] . Therefore, we 
have to add the proper momentum diffusion for atoms in 
a three-dimensional (3D) standing wave as, e.g., derived 
in Ref. .59 and given by 7 ff = {k L l at ) 2 lsc V /h\S\, where 
7 S0 is the natural linewidth of the transition. 

Taking all of these noise sources into account, we can 
now give the full QLEs describing the effective dynamics 
of the system: 



a m (t) 



a at (t) 



+ 1 



diff 



-F m {t) 



(31) 



:[fflat, H c ff] 



Taf 



Fat(t). (32) 



Here, the operator models thermal mechanical noise 
and is characterized by '(*')] = 5(t — t') and 

(£ f (W)) = ^m<5(t - t 1 ) 54J. In analogy to F m (t), 
the Hermitian noise operator F at (t) accounts for atomic 
diffusion and satisfies (F at (t)F at (t')) = — 0- 

Finally, we comment on photon loss on the optical path 
between membrane and atoms. Such loss could be mod- 
eled by placing a beam-splitter of suitable reflectivity be- 
tween the two systems, such that the light fields 6 W are 
coupled to additional vacuum modes. Integrating out the 
light fields as above would then yield a reduced coher- 
ent coupling, analogous to the situation treated in Refs. 
[481 151) . There, a similar setup without cavity was dis- 
cussed, and photon loss through the membrane into the 
— z direction was included (cf. Fig. [I]). 

IV. DISCUSSION AND APPLICATIONS 

In the previous sections, we have reduced the descrip- 



tion of our setup to the effective QLEs (31) and (32), 
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which feature a coherent coupling between atoms and 



membrane as described by the effective Hamiltonian ( 26 1 . 



However, they also contain a number of thermal and 
light-induced noise sources, whose impact we discuss in 
the following for two different scenarios. Our final goal 
is to establish interesting regimes of operation for the 
proposed setup. 



A. Coherent Dynamics 

As a first scenario we consider the observation of the 
coherent dynamics induced by the interaction term in the 
Hamiltonian ( 26 ) . A signature of the coupling is the nor- 
which is visible in the fluctuation 



mal mode splitting • 
spectra of the oscillators, as soon as it exceeds their am- 
plitude decay rates (see Ref. [5U] for an analogous discus- 
sion in standard optomechanics) . Note that normal mode 
splitting also occurs for classical oscillators and that it is 
visible also at high temperatures [50]. This effect could 
be observed with an additional weak probe laser coupled 
to the cavity that reads out the membrane vibrations. 
Analyzing the frequency spectrum of the reflected light 
would reveal the normal modes. In the time domain, the 
same effect gives rise to an oscillatory exchange of excita- 
tions between the membrane and the atomic vibrations. 
These oscillations could be observed with the probe laser 
after exciting the membrane with a piezo at frequency 

Ultimately, one would like to observe the effects of the 
coherent coupling on a single quantum level. In this con- 
text, the regime g <C w m w w a t appears naturally, where 



the effective interaction in Eq. ( 26 ) can be approximated 



by a beam-splitter Hamiltonian according to 



H, 



off 



^ m a m a m + hw at al t a at - ^r(a m a at + a a t a m) 



(33) 



which enables the coherent transfer of single excitations 
between atoms and membrane. Since all decoherence 
mechanisms are detrimental in this case, the coupling 
has to exceed all of the associated rates as expressed by 
the strong coupling conditions 



■J (at ' Im ' Im ■ 



(34) 



where 7^ = 7m N m . In order to identify a regime where 
all of these inequalities can be fulfilled, we plot in Fig. [3] 
the ratios 7^/5, 7 m lff /g, and j^/g as a function of fi- 
nesse, while fixing all other system parameters to the 
values of Tab. [TJ Note that g oc J-, such that the rel- 
ative importance 7f t lff jg of the atomic decoherence de- 
creases as 1/J-, although 7 d t lff is constant. In contrast, 
the mechanical radiation pressure noise grows quadrati- 
cally with J 7 , such that 7 m lff /<? oc T . Finally, 7^/3 scales 
as for small finesses and then saturates due to the 
laser-induced heating of the membrane for large finesses. 
The plot shows that the atomic heating is irrelevant for 




10 
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3. Coherent dynamics: We plot the ratios 'ytt 1 /g, 
lm"/9i an d "lra/9 as a function of finesse to illustrate the 
strong coupling conditions in Eq. (341. One clearly observes 
a trade-off between mechanical heating (7*) and radiation 
pressure noise (7^j a ). Parameters other than T are taken 
from Tab. HI 



the chosen parameters (which is due to the chosen large 
detuning 8), whereas there is a clear trade-off between 
mechanical heating and radiation pressure noise, with the 
latter dominating at high finesses. To illustrate this, the 
plot also shows the sum T/g of the above ratios, where 
r = 7ft ff + 7m ff + 7m ■ While T/g < 10 over the whole 
plot range, one clearly sees that there is an optimal value 
of the finesse at T ~ 300, where T/g ss 1. Interestingly, 
the condition g 3> 7 m lff implies an asymmetry in the cou- 
pling constants according to g at 3> g m (cf. Eqs. (27) and 



The above estimates show that observing normal mode 
splitting, as well as quantum coherent dynamics, is pos- 
sible for present-day parameters. Going further, the in- 
teraction in Eq. ( 33 ) can be the basis for a quantum 



state transfer protocol between the two systems (see, e.g., 
Ref. [ST]), which could become feasible for a more opti- 
mized setup. 



B. Sympathetic Cooling 

The versatility of the atoms allows one to switch on 
dissipation, e.g., by laser cooling their motion [S"Tl 152). 
Together with the resonant phonon exchange between 
mechanics and atoms contained in the Hamiltonian in 



Eq. ( 26 ) , this gives rise to a sympathetic cooling effect 



for the membrane in analogy to optomechanical cavity 
cooling. To describe this effect, we model the laser cool- 
ing by introducing an amplitude decay rate 7^° ol /2 for 
the atomic center-of-mass mode, such that the QLE (32) 
is replaced by 



Oat(*) 



. .cool 



[a at ,H eS ] - a at (f) 



^^^atW + V^^oolW. 



(35) 



Here, F coo i is a bosonic delta-correlated vacuum noise 
operator with (F cool (t)Fl ol (t')) = 6(t - f ) [Mj- 

To gain insight into the cooling dynamics, we first con- 
sider the adiabatic limit g <C 7at ol i where the atomic 
decay is so fast that the atoms adiabatically follow the 
membrane. In this case, we may eliminate the atoms and 
obtain an effective QLE for the mechanical resonator, 
from which we calculate the EOM for its mean occupa- 
tion n = (at, a™): 



n(t) = - 



cool 



[n(t) - n s 



(36) 



Here, we have introduced the sympathetic cooling rate 

.2 1 



-^cool 



v cool 



(37) 



as well as the steady state occupation of the mechanical 
oscillator 



7miVm + 7r/2 , 

7m T Tcool 

n ss i + n ss o + n ss ,3 



cool 
/at 



7r 



(38) 



The 



represent the contributions due to mechanical 
heating, 'rotating' terms in the coupling (a m a at , a^a^), 
and atomic heating, respectively, and we display them 
in Fig.E a) as a function of finesse for the parameters 
of Tab.jl It is clear that the total occupation number is 
dominated by the mechanical noise n SSl i, which we thus 
discuss in detail in the following. Note that we have 
Tcooi oc F 2 , while 7 m lff oc J 72 , and the laser-induced con- 
scales as J ' . For j m <C r cool we thus 



tribution to N n 
obtain 



27^°' 



Nuj^ t mc 2 



> 2 



(39) 
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FIG. 4. Sympathetic cooling of the membrane: Mechanical 
steady state occupation number, (a) Approximate total occu- 
pation number n ss with its individual contributions in the adi- 
abatic limit (Eq. ( [38] )) as a function of J" and for 7^°°' = 240 
kHz. For comparison, the exact mechanical occupation n™ act , 



which results from solving the QLEs ( 31 \ and ( 35 1, is depicted 
as a black line, (b) Exact steady state occupation number 
n,ss act as a function of T and 7a° o1 . The top axis displays 
the corresponding coherent coupling g and the contours indi- 
cate the associated occupation numbers. For both plots the 
parameters from Table [I] have been used. 



where a — 7r/cBMc 2 o 2 n 7 m /2fiwL|r m | 2 i4Tth and b = 
7at ol M7r 2 7 m /4w 2 t m|t m | 2 A r . Here, the three contribu- 
tions are due to light-induced diffusion, laser-induced 
heating, and non-zero support temperature, respectively. 
With the prefactors in the above expression fixed by the 
considerations of Sec.[Vj we conclude that within the adi- 
abatic limit it is advantageous to work with high finesses. 

While the above argument yields simple analytic ex- 
pressions that can be used to guide the search for op- 
timal parameters, it is only valid as long as g < 7a£ o1 . 
However, to optimize the parameters we need a solution 
also valid for larger coupling g, which can be obtained 
by solving the QLEs (31) and (35 1 exactly. For small 



finesses T < 500 and fixed laser cooling rate the result- 
ing mechanical occupation number n™ act is displayed as 
a black line in Fig.|4ja), showing that the deviations of 
the two approaches in this parameter range are not too 
large. For higher finesses, a more complete picture can 
be gained from Fig.Qb), where we show the full solution 
as a function of finesse (coupling) and cooling rate. 



Clearly, increasing the finesse beyond a certain opti- 
mal value does not yield lower occupation numbers, in 
contrast to the tendency derived above. However, the 
optimal value of the finesse roughly occurs for g ss 7at°\ 
that is on the border of the adiabatic regime. In particu- 
lar, we find n c s * act < 1 for T m 450 and 7^ Ql 220 kHz. 
We thus conclude that ground state cooling of the mem- 
brane is possible for realistic experimental parameters. 

Besides the formal similarity of our cooling scheme to 
optomechanical cavity cooling, there are two important 
differences. First, the atomic dissipation is tunable, i.e., 
after cooling the membrane close to its ground state one 
can switch off the dissipation and then observe coher- 
ent dynamics between the two systems (Sec. IV A I. This 



is not the case in cavity optomechanics, where the cav- 
ity decay is fixed. Second, our scheme does not require 
'resolved sideband' conditions for the cavity. On the con- 
trary, we want to operate in a regime where k S> w m , cj at , 
which greatly relaxes the experimental requirements on 
the cavity. In this sense, the cooling scheme described 
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here is complementary to standard optomechanical cav- 
ity cooling. 



V. EXPERIMENTAL REALIZATION 

In order to experimentally realize this proposal, we in- 
tend to use a micromechanical membrane in a Fabry- 
Perot cavity in the so-called 'membrane-in-the-middle' 
geometry [16j [55] . We will pre-cool the membrane-cavity 
setup in a cryostat where high mechanical quality factors 
up to 10 7 [B2] as well as low absorption of 10~ 6 in the 
near infrared |63j have been observed. In particular, we 
aim to use a stoichiometric SiN membrane of dimension 
1 mm 2 with a thickness of d m = 50 nm. Such a mem- 
brane has an effective mass of M = 3.6 x 10 -11 kg and a 
fundamental mechanical mode near u> m = 2ir x 400 kHz. 

On the atomic side, we consider using an ultracold 
gas of 87 Rb atoms. Because the coupling g scales as 
V~N, obtaining large atom numbers is crucial to reaching 
the quantum coherent regime and achieving ground state 
sympathetic cooling. Atom numbers as large as 3 x10 s 
have been prepared in the ground state of a large volume 
3D lattice using Raman sideband cooling |52) . We envi- 
sion using a far detuned two-dimensional (2D) lattice to 
provide confinement in the transverse direction in addi- 
tion to the longitudinal lattice that provides the coupling 
to the membrane. The resulting 3D lattice isolates the 
atoms from each other, mitigating trap loss due to light- 
assisted collisions [55]. Choosing blue detuning for the 
coupling lattice reduces spontaneous photon scattering, 
allowing for reduced laser detunings. This enables lower 
input powers to achieve the same trap frequency so that 
larger values of the cavity finesse can be obtained with- 
out causing too much absorption-induced heating of the 
membrane. 

Using the expressions for couplings and decoherence 
processes (including membrane heating) , one can find ex- 
perimental parameters that optimize the ratio of total de- 
coherence to coherent coupling or the ground state mem- 
brane occupation. In doing so, a clear trend towards min- 
imizing the atomic detuning emerges. However, such a 
prediction neglects the multi-level structure of the atom. 
For a sufficiently large lattice detuning and linear polar- 
ization, all ground-state hyperfine mp levels experience 
the same trapping potential 53 j. At smaller detunings, 
the potential becomes state-dependent and optical pump- 
ing between different mp levels induced by scattering of 
photons from the lattice light will reduce the coupling 
to the membrane. These effects can be reduced by tai- 
loring the lattice polarization or choosing a laser cooling 
technique that spin-polarizes the atoms |52j . In order to 
limit such effects, we consider 8/2n = 1 GHz as a com- 
promise. At the position of the atoms, we set the waist of 
the beam to 350 /im which corresponds to A = 9.6 x 10 -8 
m 2 . Fixing the value of 8 and A fixes the value of the 
input laser power so as to match the trap frequency to 
the mechanical frequency of the membrane. This leaves 



J- as the remaining experimental parameter, which can 



be chosen by following the discussion in Sees. |IV A| and 
|IVB| As an example, we provide an optimized param- 
eter set for u m /2ir — 400 kHz and 8/2ir = 1 GHz in 
Table |lj Values for N, Q m , and have been measured 
previously [HI [52l [621 [63] . The value of K t h is more un- 
certain as it is related to the heat conductivity of the 
membrane, which depends strongly on temperature, ma- 
terial composition and other factors. We have chosen a 
value consistent with available literature (see Appendix 

The inhomogeneous intensity profile of the lattice laser 
beams will lead to some spread Au at of atomic vibra- 
tional frequencies across the lattice, leading to dephasing 
of the atomic center-of-mass vibrations. This effect can 
be neglected as long as Aw at <C 7a° ol while laser cooling 
is applied, or Aw a t <C g in the case of coherent dynamics. 
For the numbers in Table [I] both conditions correspond 

tO Aa; a t/w at « 0.1. 
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214 kHz 


iT 


60 kHz 


<f 


8 kHz 


■7 th 


73 kHz 


5T 


4 K 





TABLE I. Optimized parameter set based on the considera- 
tions in Sec. lIVI and Sec.[V] The two left-most columns show 
mechanical and atomic parameters, respectively, while the 
column on the right displays the resulting coupling and deco- 
herence rates. 



VI. CONCLUSIONS AND OUTLOOK 

In this work, we have discussed a hybrid quantum sys- 
tem, where a micromechanical membrane in a cavity is 
coupled to the center-of-mass motion of a distant atomic 
ensemble. Placing a cavity around the membrane en- 
hances the coupling compared to the setup demonstrated 
in Ref. [51] . while preserving the modularity. Apart 
from observing coherent dynamics, the setup also enables 
bringing a cryogenically pre-cooled membrane close to its 
ground state by making use of the tunable atomic dissi- 
pation. 

In addition to cooling and manipulating micromechan- 
ical membranes, the proposed setup could also be used 
for manipulating other objects, such as levitated dielec- 
tric nanospheres [54, 65j, or even molecules. Finally, cou- 
pling to the internal atomic degrees of freedom [46l [61] 
instead of the motional ones could open additional pos- 
sibilities for manipulation and cooling, and might also 
enable even higher couplings. 



10 



ACKNOWLEDGMENTS 



Eq. ( 23 ) we obtain 



We thank A. Faber for helpful discussions. This work 
was supported by the Swiss National Science Foundation 
through NCCR QSIT and NCCR Nanoscale Science, the 
Austrian Science Fund (FWF) through SFB FOQUS, the 
DARPA OLE program, and by the EU networks AQUTE 
and MALICIA. KH acknowledges support through the 
centre of excellence QUEST, and MTR is supported by 
an EU Marie Curie International Incoming Fellowship. 



Appendix A: Treatment of Suspended End Mirror 

For the sake of completeness we also treat the standard 
optomechanical setup consisting of a Fabry-Perot cavity 
with a movable end mirror. The refractive index profile 
modeling this setup is given by Eq. (|6| with n m = 0, such 
that the mode functions read 



, , , 1 toe** 



- ik. 



) if0<z<£ 
if z > L + d 



(Al) 



Here, the cavity response is given by T' u = T L 



uj I n m — 



ie~ lKa /{l — xe 2lkL ), which for |r| — > 1 has exactly the 
Lorentzian resonances given in Eq. (12). For small dis- 
placements of the end mirror, the coupling of the latter 
to the light field can be modeled by the Hamiltonian 



H' 



A 

fk) 



(A2) 



as obtained by evaluating Maxwell's stress tensor for 
an ideal mirror, which is assumed here for simplic- 
ity. Linearizing the interaction Hamiltonian given in 
Eq. (A2) analogously to Sec.IIBl provides us with the 



same expression as in the case of the membrane given 
in Eq. (17), but with a modified coupling constant g' m = 
2afc L ; m J"/7r 3 / 2 . 



Appendix B: Elimination of the Light Field 

We briefly describe the elimination of the light field 
to obtain effective EOMs for the atomic and mechan- 
ical degrees of freedom. To this end we start from 
the Heisenberg EOMs for the quadrature operators de- 
fined by Xj = (a,j + aj)/y/2 and pj = i(aj — aj)/y/2, 
where j = {m, at} as well as x u = + b u )/V2 and 
Plj = i(bui — ^w)/v / 2- From the linearized Hamiltonian in 



X m 


= W m Pm 






(Bl) 


Pm 




f duj 




(B2) 


*^at 




25at / r— Xu, sin 

J V 27T 




(B3) 


Pat 


= — <^at#at 


[A w r + Aip u ] 


(B4) 


'f-'OJ 






(B5) 






2^m*Em 










- 2g at a; a t sin[A w r + 


A<M 


(B6) 



where again r = z/c is the propagation time between 
the two systems, A u = lu — cjl, and A(p u = <p u — (p Uh . 



Solving the Heisenberg Eqs. (B5) and (B6) for the field 
quadratures x u yields 

x u (t) = cos [A u (t - t )]x u (t Q ) + sin [A u (t - t )]p u (t ) 

-2 [ dssin[A u (t- s)] (B7) 

X 9mXm ( s ) + 5at^at(s) sin[A w r + AspJ\ , 

where the first line contains the initial conditions at some 
time t$. Under the 'bad cavity' assumption made in the 
main text, the phase difference Atp^, varies slowly over 
the frequency range of interest. Therefore, we may ex- 
pand it to first order around the laser frequency accord- 
ing to Aip URiUlj ~ A u /k, where we have assumed that 
the laser is resonant with a cavity resonance and that 
|r| — >• 1. Note that these assumptions are equivalent to 



the ones we made to obtain the Lorentzian in Eq. (12) 
In the following we thus approximate 



A<^ 



A, 



(B8) 



with the modified retardation time f = t+1/k consisting 
of the propagation time between the two systems r as well 
as the cavity lifetime 1/k. 

To proceed, we eliminate the field quadratures x u by 
inserting the solution in Eq. (B7) into Eq. (B2), which 



yields the following EOM for the mechanics 
p m = -u) m x m (t) + \ff^F m (t) 



(B9) 



2ffmffat / dsx at {s) [6g(t + T - S) - 5e(t - T - s)] 



Here, F m (t) is a noise term defined below in Eq. (B12), 
and 5e(t — s) = J^^ g |^ e - lA "('- s ) i s a representation 
of the (5-function peaked on a time-scale ~ 1/9. The 
diffusion rate of the mechanical oscillator is defined as 
7^ = 2g m . The EOM for the atomic quadrature oper- 
ator p at are obtained analogously by inserting Eq. ( |B7[ ) 
into Eq. (lB4|. 
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Summarizing these results leads to the following QLEs 
with retarded interactions: 



Pm (t) = -u m x m (t) + gx at (t -?) + fi**F m (t) (BIO) 
Pat(t) = -w at x at (t) + gx m (t - f ) + V2g at F at (t), (Bll) 

where g = 2g at gm denotes the effective coupling between 
the mechanical oscillator and the center-of-mass mode 
of the atoms, and the atomic noise term F at (t) is given 
below in Eq. (B13 1. In the above equations we can neglect 



the retardation in the coupling terms since the system 
operators evolve with rates u at , u m and g much slower 
than f _1 as indicated by the hierarchy of timescales given 
in Eq. ( 24 ) . Up to the atomic noise discussed below we 



thus recover the EOMs given in ( 25 ) of the main text 



The noise operators appearing in the above equations 
are Hermitian and given by 

F m (t) = -y/2 J -|L { cos [A u (i - t )]x u {t ) 

+ sin[A LJ (t-io)]Pc(<o)} (B12) 

F at {t)=-^j -jL sin[A w f]{cos[A w (t-i )K(to) 
+ sm[A u ,(t-to)]pu,(to)}- (B13) 



To obtain the corresponding correlation functions, we as- 
sume that the field modes are initially in the vacuum 
state and therefore obey (&u;(io)fyL'(io)) = S(u — oj'). 

For the correlation function of the mechanical noise 
term we obtain (F m (t)F m (t')) = S g (t ~ if). Note that in 
the absence of the atoms our setup is exactly the one con- 
sidered in standard cavity optomechanics, and we should 
thus find similar decoherence mechanisms. In order to 
compare the above diffusion rate 7 n i 1 lfI to the optomechan- 



diff 



(go<=»c) 



ical cooling literature we can express it as 7J 
where we used that the intra-cavity amplitude is given by 
a c = a/y/2n and that the single photon optomechanical 
coupling can be expressed as go = fj^m = 2|r m |w c / m /L 
[16] . When taking into account that we assumed the 'bad 
cavity' regime where n 3> w m and that we drive the cavity 
on resonance, we obtain the same scaling of the diffusion 
rate as presented in [66 , which was recently observed in 
an experiment |58j . 

The correlation function for the atomic noise term 
reads {F at {t)F at {t')) = -\[S e (t+2f-t')+Sg{t-2f-t')- 
25g(t — t')]. On the timescale of the atomic evolu- 
tion the three 5-kicks in the correlation function average 
to zero due to our assumption of small retardation times 
(w at f <C 1, see Eq. (24)). For this reason the atomic 
noise operator F at (t) does not occur in the QLEs in the 
main text in Eqs. ( [25] ) However, this is an artifact of the 
ID treatment, since atoms in a standing wave are well- 
known to undergo momentum diffusion |59j . This effect 
is introduced in Sec. lIII Bl 



Appendix C: Membrane Heating Due to Laser 
Absorption 

In this section membrane heating due to laser absorp- 
tion is explicitly modeled. For a circular membrane of 
diameter I, a simple analytical solution can be obtained 
and the case of the square membrane can then be ap- 
proximately described by including a 'geometric prefac- 
tor' of order unity. The results of the analytical model 
are checked against a finite element method (FEM) sim- 
ulation. 

Heat transport inside the SiN membrane is governed 
by the heat equation for the temperature T: 



dT 

PCp-fo = K thAT + Q th . 



(CI) 



where Kth is the thermal conductivity, p the mass density, 
c p the specific heat capacity of the membrane material, 
A the Laplace operator, and Qth is a source term describ- 
ing the power dissipated per unit volume. For a 2D tem- 
perature distribution inside a thin membrane (thickness 
dm *C I) in steady-state (dT/dt = 0), the heat equation 
reads 



AT=--(r — 
r dr V dr 



1 d fdT 



r 2 dip \ dcj) 



= -Qth/K th . (C2) 



To obtain a simple analytic solution, the case of a circular 
membrane placed in the center of the intra-cavity elec- 
tric field mode is considered. For simplicity, we assume 
that the power P a b s is absorbed within a circle of area 



•Art 



TTW„ 



where u> m is the waist (e radius) of the 



Gaussian field mode at the position of the membrane. 
Thus, Qth = -Pabs/(-4mdm)- The absorbed power for a 
membrane sitting on a slope of the intra-cavity intensity 



standing wave is given by P abs = a^ n 2P circ = a, 2 n ^P 
with the power absorption coefficient a^ n = 1 — — t^. 
Here, P c i rc is the intra-cavity circulating running wave 
power defined by P c i rc = 2 ^-P and P is the power of the 
incident laser. The peak laser power for a standing wave 
inside the cavity is given by P s t = 4P c i rc , such that we 
obtain half peak power at the position of the membrane. 
Making use of the azimuthal symmetry (dT/d<fr = 0), the 
equation simplifies to 



1 d ( dT 



r dr V dr 



-Qth/^th for < r < w m , 
for w m < r < 1/2. 



(C3) 



We solve this equation subject to the boundary condition 
that the frame is at constant temperature, T(r = 1/2) = 
Tq, yielding the temperature distribution 



T(r) 



'T + <5T- ¥^r 2 



for < r < 



^o + %^ln(i) for w m <r< 1/2, 
where 

5T = T{r = 0) - T 



2«th 



In 



(C4) 



(C5) 
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is the temperature increase of the membrane center com- 
pared to the frame. The average membrane temperature, 
obtained by integrating T(r) over the membrane, is 



To 



4« t h 



1 



(C6) 



The thermal link K t h connecting the membrane center 
to the frame is defined as K t \ x = P a ^ s /5T. Using P a b s = 
QthTU^dm, we find 



K, 



th 



-^abs 

~5T 



2TTK t hd n 



In 



(ihL) 



(C7) 



The values for ST and -Kth calculated above for a circu- 
lar membrane are confirmed by a FEM simulation. Run- 
ning the FEM simulation instead for a square membrane 
of side length I, keeping all other parameters fixed, we 
obtain similar results but with a geometric prefactor f g 
so that 



square — fg 8-11(1 ^th, square — ^th/ ' fg 



(C8) 



are the corresponding quantities for the center of the 
square membrane. We find that f g decreases from f g — 
1.075 for w m /l = 0.3 to f g = 1.017 for w m /l = 0.01. This 
shows that the model of a circular membrane gives an 
excellent approximation for the square membrane over 
the relevant parameter range, well within the errors of 
other parameters (such as K t h)- Note that in general 
w m /l < 0.3 to avoid clipping of the beam at the mem- 
brane frame. 



As shown in Eq. (C7|, the thermal link K t h depends on 



the membrane geometry, the beam waist, and the ther- 
mal conductivity Kth- Because Kth depends strongly on 
temperature, material composition, and other factors its 
value is not precisely known. In [67] , a room-temperature 
value of K t h = 25 — 36 WK _1 m _1 is given for stoichiomet- 
ric SiN grown by chemical vapor deposition. We extrap- 
olate this value to low temperatures using a temperature 
dependence as in |68j . 
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